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We have theoretically studied successive transitions in vanadium spinel oxides with (£29) 2 
electron configuration. These compounds show a structural transition at ~ 50K and an an- 
tiferromagnetic transition at ~ 40K. Since threefold ti g orbitals of vanadium cations are 
occupied partially and vanadiums constitute a geometrically-frustrated pyrochlore lattice, 
the system provides a particular example to investigate the interplay among spin, orbital and 
lattice degrees of freedom on frustrated lattice. We examine the models with the Jahn-Teller 
coupling and/or the spin-orbital superexchange interaction, and conclude that keen compe- 
tition between these two contributions explains the thermodynamics of vanadium spinels. 
Effects of quantum fluctuations as well as relativistic spin-orbit coupling are also discussed. 

§1. Introduction 

In this paper, we review our recent studies on phase transitions in vanadium 
spinel compounds, AV2O4 (74=Zn, Mg or Cd). 1 ^ 5 ) This family of materials is an- 
tiferromagnetic insulator with spin S = 1 on pyrochlore sublattice in the spinel 
structure (Fig. P). Recent experiments showed the presence of two successive phase 
transitions at low temperatures around 50K. 6 ) As temperature is lowered, a struc- 
tural transition occurs first, then a magnetic transition follows. As we will discuss 
later, simple spin exchange mechanism does not explain the lowest temperature 
phase, and we need to take account of two important factors in order to explain 
basic character of these transitions. One is geometrical frustration and the other is 
strong correlation effects, which lead to keen competition of various degrees of free- 
dom. As widely known, the geometrical frustration often leads to huge degeneracy 
of the ground-state manifold, and suppresses a long-range ordering. On the other 
hand, the strong electron correlation opens several channels in the low energy sector, 
such as spin and orbital degrees of freedom. A key question is how these channels af- 
fect the manifold — do they result in some nontrivial phenomena through enhanced 
fluctuations, or lift the degeneracy in some nontrivial manner? We would like to 
clarify this intriguing problem by studying effects of both frustration and correlation 
in the vanadium spinels. 

Vanadium spinels are a 3ti-electron system with orbital degrees of freedom. In 
the last decade, various transition-metal materials have been intensively studied 
to investigate the effects of orbital degrees of freedom. 7 ) In particular, manganites 
showing colossal magnetoresistance (CMR) are the most typical example. Vanadium 
spinels have different features from these materials. A distinct difference is the lattice 
structure. In manganites showing CMR, the lattice structure is perovskite; atomic 
bonds are connected with right angle and the network of magnetic ions is close to a 
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Fig. 1. Cubic unit cell of the pyrochlore lattice. In vanadium spinels, vanadium cations constitute 
this pyrochlore structure. 

simple cubic lattice. In spinels, atomic bonds are connected to constitute a lattice 
in which elementary plaquettes are triangles; the network of magnetic cations is the 
pyrochlore lattice shown in Fig. ^ Another important difference is that the vana- 
dates are a t2g~electron system, whereas the manganites have conduction electrons 
in e g multiplet. Jahn-Teller distortions are generally smaller in ti g systems, and the 
crystal-field splitting of the electronic levels is less prominent, partly because ti g or- 
bitals point away from surrounding anions. Therefore, it is important to take account 
of orbital degrees of freedom in ti g systems like vanadium spinels. We also note that 
tig electrons may have a nonzero orbital angular momentum, whereas orbital angular 
momenta are quenched in e g systems. This also leads to several effects present only 
in tig systems. For example, relativistic spin-orbit coupling is sometimes important 
to take into account, but this is not a central issue in this paper. 

From the viewpoints of orbital physics and geometrical frustration, spinels with 
transition-metal ions are very interesting system. In this paper, we shall discuss 
about two phase transitions in vanadium spinels with divalent A-site cations such as 
Zn, Mg or Cd. In these materials, first of all, interesting phenomena such as phase 
transitions occur at much lower temperatures compared to a typical energy scale in 
the system. That is, the characteristic energy scale, which is given by, for example, 
the Curie-Weiss temperature estimated from the temperature dependence of the 
magnetic susceptibility, is ~ 1000K, 8 ) nevertheless the system remains paramagnetic 
and does not show any transition down to ~ 50K. This suppression of relevant energy 
scale is characteristic to frustrated systems. At low temperatures, in the end, the 
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Fig. 2. Magnetic ordering structure in the lowest-temperature phase of vanadium spinels deter- 
mined by neutron scattering experiment. 9 - ) 

compounds show two successive transitions. 6 ) One is a structural transition at ~ 50K 
and the lattice changes its symmetry. While the lattice symmetry is cubic above 
the transition temperature, it changes to tetragonal below the temperature, and the 
system is compressed along one of the principal axes, the c axis. The other transition 
takes place at a lower temperature ~ 40K, and ascribed to an antiferromagnetic 
ordering. The magnetic structure in the lowest-temperature phase was determined 
by neutron scattering experiment as shown in Fig. |21 9 ) It consists of up-down-up- 
down- staggered antiferromagnetic pattern along the chains in the xy planes [(110) 
and (110) directions] and up-up-down-down- pattern with period four along the 
chains in the yz and zx planes [(101), (101), (Oil) and (Oil) directions]. Thus 
the fundamental questions are (i) what the driving mechanism of these successive 
transitions is and (ii) how the complicated magnetic ordering is stabilized. We will 
answer these questions theoretically by investigating keen competitions among spin, 
orbital and lattice degrees of freedom in this geometrically-frustrated system. 

This paper is organized as follows. In Sec. [21 we examine a model with the Jahn- 
Teller coupling alone. We also examine a model with spin-orbital superexchange 
interaction alone in Sec. El In these two sections, we show that both two models 
fail to explain the experimental results. In Sec. |11 we take into account both two 
terms in the model and study keen competition between them. We show, by using 
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Fig. 3. Crystal-field splitting of d levels. 



mean-field analysis and Monte Carlo simulation, that the model successfully explains 
two successive transitions of vanadium spinels. Effects of fluctuations, both thermal 
and quantum, are also discussed. In Sec. 03 we make some remarks on effects of 
relativistic spin-orbital coupling, direction of magnetic moment, and a possibility of 
Haldane gap state in this system. Section H3 is devoted to summary. 

§2. Model with Jahn Teller coupling 

Let us start from a simple model with the Jahn-Teller coupling alone, which 
describes a coupling of orbital degrees of freedom with lattice distortions. In the 
octahedral position surrounded by a symmetric Oq octahedron, the cubic crystal 
field splits the fivefold d levels of vanadium cation into higher twofold e g levels and 
lower threefold tig levels (left two columns in Fig. |3J). Jahn-Teller distortion is a 
distortion of Oq octahedron, which gains energy by lifting the degeneracy further to 
achieve a closed-shell configuration of electronic state. In other words, the Jahn- 
Teller effect is active only for the case of an open-shell configuration, and hence, 
depends on the number of d electrons. 

In vanadium spinels with divalent ^4-site cation ^4V2C>4, each vanadium has 
valence V 3+ , and this corresponds to (i2g) 2 configuration. Hence, threefold tig levels 
are partially filled (second column in Fig. |3J) so that the Jahn-Teller effect is active. 
We start from an ionic model describing independent VOq octahedra, and later 
we will introduce intersite correlations. If we consider the Jahn-Teller distortion of 
tetragonal symmetry, which is observed in experiments, the ionic model can be given 
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in the form 

#j T ic = E ^fatt + n l2 - 2n i3 ) + \Y,Ql I 2 " 1 ) 

i i 

where 7 is the electron-phonon coupling constant and Qi denotes the amplitude of 
local lattice distortion at site i. Here, we consider the tetragonal distortion along one 
of the principal axes, the z direction, and choose the sign of Qi such that it is positive 
for a compression of the octahedron in the z direction. The orbital indices a = 1, 2 
and 3 of the density operator ni a denotes d yz ,d zx and d xy , respectively. Assuming 
high-spin states by strong Hund's-rule coupling, the model is defined for spinless 
fermions. The second term describes the elastic energy (Q is normalized to absorb 
the elastic constant); here we use the adiabatic approximation and phonons are 
regarded as classical objects. The model (|2- X|> is considered with the local constraint 
on the density as 

]T(n ia > = 2. (2-2) 

a 

The ground state of the ionic model (|2T[) is obtained by minimizing the energy 
under the constraint (|2-2|) . The ground-state energy per site is given by 

£i° T nic = -2 7 2 for Qi = -2 7 , (2-3) 

and for the orbital states 

(nil) = (n i2 ) = 1 and (n i3 ) = 0. (2-4) 

This situation corresponds to the elongation of VOg octahedra as shown in the third 
column in Fig. |3J The doublet of d yz and d zx levels is lowered and occupied by two 
electrons to form the closed-shell configuration. 

Experimentally, lattice distortion occurs in the opposite way at low tempera- 
tures, that is, a uniform compression of octahedra. If we restrict to the case of 
compression Q > 0, the energy minimum is found as 

£j°T nic = -y for Q i = 7 , (2-5) 

and for the orbital states 

(^1 + n i2 ) = 1 and (n i3 ) = 1, (2-6) 

which are shown in the last column in Fig. |H1 

In the real lattice structure, the Jahn-Teller effect should have some cooperative 
aspect — V06 octahedra share oxygen ions between neighbors, which leads to inter- 
site correlations of the Jahn-Teller distortion. Since spinel lattice structure consists 
of three-dimensional network of edge-sharing octahedra as shown in Fig. 21 a distor- 
tion of VC>6 octahedron tends to distort neighboring octahedra in the same manner. 
In other words, a ferro-type correlation is expected for the tetragonal distortions be- 
tween neighboring octahedra. This tendency is mimicked by the ferro-type intersite 
coupling of the distortion as 

^-A^^, (2-7) 
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Fig. 4. Edge-sharing network of VOe octahedra in the spinel structure. Vanadium cations locate 
in the center of each octahedron and constitute the pyrochlore lattice in Fig. Q 



where A > is the coupling constant and the summation is taken over the nearest- 
neighbor sites on the pyrochlore lattice which consists of vanadium cations. 

If we omit farther-neighbor correlations beyond the nearest neighbors, the total 
Hamiltonian for the Jahn-Teller effect is given by 

U zjionic | j_rnn /r> q\ 

We can minimize the total energy with respect to Qi by assuming uniform orbital 
states at all sites. This corresponds to find a mean-field solution for cooperative 
Jahn-Teller distortion of ferro-type, which is naturally expected for the ferro-type 
intersite coupling in Eq. (|2-7|) . For the orbital state (|2-4|) . or equivalently for an 
elongation of all V06 octahedra, we obtain the minimum of energy per site as 

Note that 0<A<l/6is required to keep the system stable. On the other hand, 
for the orbital state (|2-6|) . corresponding to a compression of octahedra, we find a 
minimum as 

^ T = "2(l36A) fOT * = r?6V (24 °» 

Since Ejt < -Ejt, the model with Jahn-Teller effect alone predicts a uniform elonga- 
tion of all V06 octahedra, which apparently disagrees with the experimental result. 
In order to stabilize the compression of octahedra as observed in experiments, we 
need to compensate the energy difference per site 

3 7 2 

A]T = £ JT _ £ JT = — — L— — (2-11) 
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by some energy gain from other mechanisms. We will see in Sec. 0] that this com- 
pensation indeed comes from the orbital superexchange interaction. 

§3. Model with superexchange interaction 

3.1. Model Hamiltonian 

Next we examine effects of spin and orbital superexchange interactions. The 
superexchange interactions describe effective intersite couplings of both spin and 
orbital degrees of freedom in low-energy physics of strongly-correlated insulators. 
They are derived from a multi-orbital Hubbard model by using the perturbation in 
the strong-correlation limit. 

We start from the Hubbard model with threefold orbital degeneracy correspond- 
ing to the t2 g manifold. The Hamiltonian reads 

#Hub = ^(*»«,J'/9 C Lt^t + H.C.) + - ^ U aP,a>P>c\ aT c\p T ,C i( 3> T >C ia > T , 

ij a/3 r i a(3,a' (3' tt' 

(3-1) 

where i,j and r, r' are site and spin indices, respectively, and a,/3 = 1 (d yz ), 2 (d zx ), 
3 (d xy ) are orbital indices. The first term of //Hub is the electron hopping, and the 
second term describes onsite Coulomb interactions, for which we use the standard 
parametrizations, 

U a /3,a'p' = U'5 a a'S/3f3' + Ja(8 a p' S/3a' + 8 a /35a'/3'), (3-2) 

U = U' + 2J R . (3-3) 

In order to describe the low-energy effective theory for the insulating state of 
vanadium spinels, we apply the second-order perturbation in the hopping term. 
The unperturbed states are atomic eigenstates with two electrons on each vana- 
dium cations in a high-spin state (S = 1). Among various hopping processes, the 
most relevant contribution is the nearest-neighbor a bond, because the spinel lattice 
structure consists of three-dimensional network of edge-sharing V06 octahedra in 
Fig. 2] 10 )' 11 ) This is the hopping in the case where one of four lobes of each t2 g 
orbital is pointing towards the other end of the bond. In the perturbation, we take 
account of this u-bond contribution only for simplicity. The superexchange Hamil- 
tonian, which is called the Kugel-Khomskii type Hamiltonian, 12 ) is obtained for the 
nearest-neighbor sites in the form 

(3-4) 

n ia(ij)) n ja(ij)] > (3'5) 

(3-6) 

where Si is the S = 1 spin operator and ni a = ^ T c\ aT Ci aT is the density operator 
for site i and orbital a. Here, a(ij) is the orbital which has a finite hopping integral 



-"SE - J <■ ""o-AF + %-F J' 

fro-AF = ( A + BS i ■ Sj) [n ia{ij) (l - n ja{ij) ) + (1 - 

ho—F = C(l — &i ' Sj) n ia(ij) n ja(ij)i 
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between the sites i and j, for instance, a(ij) = 3 (d xy ) for i and j sites in the same 
xy plane. Parameters in Eqs. (|3-4j) - (|3-6|) are given by coupling constants in Eq. (|3- lft 

as 

J= { -jj L , (3-7) 

A=*-±, 5 = -5-, C=^L, (3-8) 
1 — 3»7 I — 377' 1 + 27/' 1 ; 

where i° n is the hopping integral of a bond, and rj = Jn/U is a small parameter of 
~ 0.1. Note that each site is subject to the local constraint, Yla=i n ia = 2. 

The superexchange model (|3-4|) consists of two contributions, h^ AF and %_ F : 
The former favors spin-ferro and orbital-antiferro configuration on ij bond while 
the latter favors spin-antiferro and orbital-ferro configuration. The competition 
between these two exclusive contributions is a general aspect of Kugel-Khomskii 
type models. Another characteristic feature is the difference in the symmetry of 
the interactions. In the spin sector, the superexchange coupling is Si ■ Sj, isotropic 
Heisenberg type and independent of the bond direction. On the other hand, the 
superexchange coupling in the orbital sector becomes anisotropic because transfer 
integrals depend on orbital indices as well as the bond direction. Particularly in the 
present model (|3-4|) . the anisotropy is conspicuous. The orbital exchange interaction 
depends only on the density operator, in other words, it is a three-state clock type 
interaction — there is no quantum fluctuation because the density operator is a 
constant of motion. Moreover, the orbital interaction depends on the bond direction 
and the orbital states in two sites in a complicated way. These features come from 
the dominant role of transfer integrals of a bond, which are orbital diagonal and 
strongly depend on the bond direction and the orbital states. 

3.2. Ground state 

Now we examine the ground state of the superexchange Hamiltonian (|3-4|) de- 
fined on the pyrochlore lattice. Here, we consider the ground state with four- 
sublattice ordering in V4 tetrahedron unit cell. Let us suppose rj = for a while. 
In this case, the energy is minimized for two different spin-orbital configurations 
in Fig. |SJ The energy per site is E^e = — 4 J, with counting the interactions with 
neighboring tetrahedra. Per tetrahedron unit cell, the configuration (a) contains four 
spin-antiferro orbital-ferro bonds (the energy is — 2J per bond), four spin-antiferro 
orbital-antiferro bonds (—«/), and four spin-ferro orbital-antiferro bonds (—</). The 
configuration (b) contains four spin-antiferro orbital-ferro bonds, two spin-antiferro 
orbital-antiferro bonds, and six spin-ferro orbital-antiferro bonds. 

The degeneracy between these two states is lifted for 77 > 0. Energies per site 
for configuration (a) and (b) are given by 



= -2J(A + C), tfg: = - 2 J ( A + \ B + C l ' : '> A) 1 



respectively. Since B is positive, the configuration (b) has a lower energy. This 
appears to be inconsistent with experiments: The spin structure in (b) is incompat- 
ible with the antiferromagnetic state observed in experiments in Fig. [5] because two 
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Fig. 5. Candidates for the ground states with four-sublattice ordering of the superexchange model 
13-411 . Spin and orbital configurations in tetrahedron unit cell are shown. Arrows denote spins, 
and labels (d yz , d xy ), etc., indicate two occupied orbitals at each site. 

neighboring spins in the lower xy plane are parallel whereas the experimental spin 
configuration is antiparallel in all xy planes. Moreover, configuration (b) noticeably 
breaks the tetragonal symmetry which is experimentally observed in structurally- 
ordered phase. Hence, this mean-field type analysis for the ground state tells that 
the model with superexchange interactions alone fails to reproduce the experimental 
results. 13 ^ 

There is, however, one caveat — the energy difference between the ground state 
(b) and the next lowest state (a) is small compared to the relevant energy scale in 
experiments. The energy difference between (a) and (b) per site is given by 



Let us roughly estimate this energy difference. As discussed in Ref. 2), estimates for 
vanadium spinels are J ~ 200K and rj ~ 0.1. Hence, the energy difference becomes 
Z\ SE ~ 30K, which is smaller than the structural transition temperature ~ 50K. This 
strongly suggests that spin and orbital configuration in the ground state is sensitive 
to competing interactions; in particular, the stable configuration may be changed by 
including the cooperative Jahn-Teller effect, which apparently plays an important 
role in the structural transition. Indeed, we will see that the energy for the state (a) 
can be lower than for (b) in some parameter regime of the Jahn-Teller coupling - 
see Sec. HJ 

3.3. Instability from high-temperature phase 

We can gain further insight into the stable spin and orbital configuration of the 
model (|3-4[) by using another mean-field type argument. We here discuss fluctuation 
effects and an instability in spin and orbital sectors when we decrease temperature 
from high-temperature phase. 1 ) 

Let us first examine the fluctuations in the spin sector. For this purpose, we 
replace the density operators in the orbital part of model ()3-4|) by their expectation 



(3-10) 



1 -3r/ 
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values in the high-temperature regime; rii a — * (rii a ) = 2/3. Then we end up with 
the reduced Hamiltonian describing the spin sector only, which reads 

H spm = Js'Y]S i - Sj, (3-11) 

up to irrelevant constants, where the effective exchange constant is given by 

J S = P(C-B). (3-12) 

Since Js is positive for small r\ ~ 0.1, the reduced model (|3-llj) is the antiferromag- 
netic Heisenberg model with 5 = 1 spins defined on the pyrochlore lattice. It is 
known that the spin correlations hardly develop in this model because of the strong 
geometrical frustration. A possibility of a valence-bond-solid type ordering with the 
aid of spin-lattice coupling has been discussed, 14 ) however it is not easy to explain 
two successive transitions in vanadium spinels by this spin Jahn-Teller scenario. 
The present consideration implies that the magnetic instability cannot occur in the 
high-temperature orbital-para phase. 

Next we consider the instability in the orbital sector. In this case, we replace 
the spin part of model ()3-4(l by their mean value in high-temperature regime, that 
is, Si — > (Si) = 0. The resultant Hamiltonian for the orbital degree of freedom reads 

-^orbital = ^ ~] n ia(ij) n ja(ij)i (3'13) 
(ij) 

up to irrelevant constants, where 

Jo = J(2A-C). (3-14) 

This coupling constant Jo is also positive so that Eq. ()3T3j) is also the antiferro-type 
model on the pyrochlore lattice as the model (|3T1|) . A crucial difference between 
the two models is the form of interactions: The orbital model (|3T3p has anisotropic 
interaction (three-state clock type) which strongly depends on the bond direction 
and the orbital states. This peculiar form of the orbital coupling lifts the degeneracy 
in this pyrochlore system, at least partially — the three-state clock model still suffers 
from a frustration on the pyrochlore lattice. The lowest energy is achieved by four 
orbital-antiferro bonds and two orbital-ferro bonds within each tetrahedron. There 
are two different types of orbital configurations with these bonds as shown in Fig. H3 
In type (a), two orbital-ferro bonds do not touch with each other, while they touch 
at one corner of the tetrahedron in type (b) . Note that these orbital configurations 
coincide with those in Fig. [5J The energy of the reduced orbital model (|3T3[) is the 
same for these two configurations. Hence, there remains degeneracy in the orbital 
sector although the degeneracy is reduced and not severe compared to that in the 
spin sector in the model (|3Tlj) . 

The degeneracy remaining in the orbital sector is easily lifted by including the 
Jahn-Teller effects. This will be discussed in Sec. EJ 
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Fig. 6. Stable orbital configurations of the effective orbital model 13-1311 . Labels (d yz , d xy ), etc., 
indicate two occupied orbitals at each site. Thick lines denote orbital-antiferro bonds. 

§4. Model with both superexchange and Jahn Teller couplings 

In this section, we examine the model with both superexchange interactions and 
Jahn-Teller couplings. In Sec. 14.11 we will show how the mean-field arguments in 
Sec. 01 are modified by including Jahn-Teller distortions in addition to the superex- 
change model. We will point out that once both superexchange and Jahn-Teller 
effects are taken into account, the degeneracy in the orbital sector is lifted and an 
orbital ordering takes place. We will also show that spin frustration partially remains 
in the presence of the orbital ordering, which is finally reduced by farther-neighbor 
exchange couplings. In Sec. 14.21 the mean-field picture is examined and confirmed 
by Monte Carlo simulation. Thermodynamic properties are compared with experi- 
mental results. In Sec. 14.31 effect of quantum fluctuations, which is neglected in both 
mean-field argument and Monte Carlo simulation, is examined by using linearized 
spin-wave theory. 

4.1. Mean-field type analysis 

We have shown in Sec. that the Kugel-Khomskii type of spin-orbital model 
has two competing stable states for V4 tetrahedron unit. For realistic values of 
parameters for the vanadium spinel, it was found that the lowest-energy state has 
a different symmetry from the experimental results and is not compatible with the 
antiferromagnetic order in experiments. This does not imply the failure of the Kugel- 
Khomskii model but that one needs to take account of couplings to lattice distortion. 

Below the structural transition temperature at ~ 50K, the lattice shows a 
tetragonal distortion, and VOg octahedra are uniformly compressed along the c axis 
(c < a). Energy level scheme in this case is shown at right most in Fig. 1, and three 
degenerate t2 g orbitals in cubic crystal field are split to higher-energy doublet and 
lower-energy singlet states in tetragonal field: 



As shown in Eq. I|2T|1 . the level splitting A = e(d yZjZX ) — e(d xy ) is proportional to 
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the tetragonal distortion, A = 3jQ with Q oc 1 — c/a. 

Let us first reanalyze stable orbital configuration in spin paramagnetic phase 
with taking account of the Jahn-Teller effects. We showed in Sec. 13.31 that two 
types of orbital configurations in Fig, H3 have the same energy of the effective orbital 
interaction (3.12). Configuration (a) is compatible with the tetragonal symmetry 
as that of the intermediate-temperature phase experimentally observed below the 
structure transition temperature, whereas configuration (b) has a lower symmetry. 
Energy gain due to the Jahn-Teller effects is largest when the lattice is compressed 
along the z axis for configuration (a) but along the x or y axis for configuration (b). 
The corresponding energies per site are respectively, 



< = --. (4-2) 
3 12 



Therefore, the tetragonal orbital configuration (a) has a lower energy and is stabilized 
by the Jahn-Teller effects, as far as a uniform mode of distortion is considered. 

Next, let us reexamine the stability at T = 0, discussed in Sec. 13.21 by including 
the effects of tetragonal lattice distortions. The energy gain of the configuration is 
largest when the lattice is compressed in the same matter as for the previous orbital 
model. The configurations shown in Fig. 2 (a) and (b) have energies, respectively, 

4e+jt = 4e + 4? = + C) - j, (4-3) 

= + < = -2J(A + \b + C) - ^, (4-4) 

where E^'^ are given in Eq. ()3-9|) . Comparing these two, one finds that the 
configuration (a) has a lower energy than (b) when 

A>4JB = -i^-. (4-5) 
l-3?7 

As was emphasized in previous sections, the vanadium spinels are the sys- 
tem in which the Jahn-Teller couplings compete with the spin-orbital interactions, 
Z\~ J x (coordination number). Recalling that rj is small, we may expect that the 
above condition is satisfied and that the configuration (a) with tetragonal symmetry 
is most stable, if spin-orbital superexchange interactions and Jahn-Teller couplings 
are both considered. In Sec. 14.21 we will demonstrate by Monte Carlo calculations 
that the tetragonal orbital order (a) actually appears at finite temperatures. 

Next, let us discuss the spin configuration in the magnetic phase. The config- 
uration in Fig. |Hf a ) imposes the spin configuration inside the tetrahedron unit cell, 
and one needs to consider correlations between different tetrahedra. The spin order 
at ~ 40K was observed by neutron scattering experiment 9 ) and the determined spin 
pattern is shown in Fig. [21 This pattern corresponds to the Q = (0,0, 2tt/c) order 
of the configuration in Fig. [5fa) (c is the length of cubic unit cell along the z axis), 
i.e., uniform packing in the xy planes and alternate stacking in the z direction with 
spin inversion from layer to layer. Antiferromagnetic spin pattern along chains in 
the xy planes are easily explained, since any pair of nearest-neighbor spins in xy 
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Fig. 7. Frustration between neighboring xy chains remaining in the mean-field level. The mean- 
field energy does not depend on the relative angle 9 between the antiferromagnetic chains. 




Fig. 8. (a) Second-neighbor exchange couplings, and (b) third-neighbor exchange couplings. 

chains should be anti-parallel. However, there are no correlations between different 
xy chains within the level of the mean-field approximation. Look the spin configura- 
tion in Fig. Efa) and rotate two spins in the lower xy plane with imposing that they 
are always anti-parallel to each other as in Fig. [Tj Such rotations do not change the 
energy in the mean-field level, and this implies no correlation between two chains. 
This is due to the frustrated lattice structure and antiferromagnetic nearest-neighbor 
correlations in each chain. 

Since interchain correlations are not determined by nearest-neighbor couplings, 
we have to consider the effects of longer-range couplings, which are much weaker in 
amplitude than the nearest-neighbor ones and are usually neglected. These longer- 
range couplings are shown in Fig. |SJ Another possible mechanism is order by dis- 
order. 15 ) However, as we will discuss later based on numerical (Sec. 14.2(1 and ana- 
lytical (Sec. 14.3(1 calculations, the mechanism of order by disorder is not sufficient 
to explain a finite static moment in the magnetic phase. As shown in Fig. once 
antiferromagnetic correlations develop in the xy chains, it is important to note that 
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Fig. 9. (a) Two sets of spin chains coupled by third-neighbor exchange couplings J3 are depicted 
in black and white colors, (b) Second-neighbor coupling J2 are strongly frustrated between 
neighboring antiferromagnetic chains in the xy planes. 



third-neighbor couplings lead to more dominant effects than second-neighbor ones. 
Third-neighbor couplings connect antiferromagnetic chains every two xy planes as 
shown in Fig. EO^a). As a consequence, the whole system becomes two interpenetrat- 
ing sets of three-dimensional spin system [black and white spins in Fig. Efa), re- 
spectively]. These two sets are connected by nearest-neighbor and second-neighbor 
couplings. However, their mean fields are canceled out and lead to no effects. The 
frustration of second-neighbor exchanges is depicted in Fig. Efb). In this case, it is 
fluctuations which determine the correlation between the two sets. We discussed the 
effects of quantum fluctuations, in particular, in Ref. 1), and it was found that the 
collinear spin order is stabilized. Monte Carlo results in the next section show that 
thermal fluctuations also stabilize the same collinear order. 

4.2. Monte Carlo study 

Now we investigate the thermodynamic properties of the model with both su- 
perexchange interaction and Jahn-Teller coupling by using Monte Carlo method. In 
Sec. 14.2.11 we briefly describe the Monte Carlo technique and the model parameters 
used in the simulation. In Sec. 14.2.21 the model with nearest-neighbor couplings is 
studied. We will show that the orbital ordering takes place at a finite temperature 
but the magnetic frustration remains down to the lowest temperature. Effects of 
third-neighbor exchange interactions will be examined in Sec. 14.2.51 As predicted in 
the mean-field argument, Monte Carlo results show that the magnetic frustration is 
reduced and antiferromagnetic order is stabilized at a finite temperature. The phase 
diagram by controlling the third-neighbor coupling will be shown in Sec. 14.2.41 

4.2.1. Method and parameters 

In order to investigate thermodynamic properties of the effective spin-orbital- 
lattice coupled model, we use Monte Carlo calculations. Quantum Monte Carlo 
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simulation is known to be difficult for frustrated models because of the negative sign 
problem. In the present study, we neglect quantum fluctuations and approximate 
the model in the classical level. This approximation retains effects of thermal fluctu- 
ations that may play dominant roles in finite-temperature transitions. The quantum 
nature originates only from spin 5 = 1 operators because the orbital interaction in 
Eq. (|3-4j) is classical, and because Jahn-Teller lattice distortions are treated as clas- 
sical variables in Eq. (|2-8|) . We approximate the spin operators by classical vectors 
with the modulus \S\ = 1 to give largest z component S z = 1 in the classical part. 
Therefore, our model for Monte Carlo simulation consists of the classical Heisen- 
berg part in the spin sector, the three-state clock part in the orbital sector, and the 
classical phonon part. 

We use a standard metropolis algorithm with local updates. We typically per- 
form 10 5 Monte Carlo samplings for measurements after 10 5 steps for equilibration. 
The system sizes in the present work are up to L = 12, where L is the linear dimen- 
sion of the system measured in the cubic units, i.e., the total number of sites N s a e is 
up to 12 3 x 16 = 27648. We set J = 1 as an energy unit and the lattice constant of 
cubic unit cell as a length unit. We use the convention of the Boltzmann constant 
£?b = 1- The following results are for rj = 0.08, 7 2 = 0.04 and A = 0.15. Readers are 
referred to Ref. 2) for more details of Monte Carlo simulation. 

4.2.2. Model with nearest-neighbor couplings only 

Let us first examine what happens in the model with the nearest-neighbor in- 
teractions only, whose Hamiltonian consists of the nearest-neighbor superexchange 
interactions (|3-4j) and the Jahn-Teller coupling (|2-8|) ; 



The mean-field type analysis in Sec. 14.11 predicts that the orbital-ordered state in 
Fig. Ela) is stabilized by a cooperative effect of orbital interaction and Jahn-Teller 
coupling, but spin frustration partially remains even in the presence of the orbital 
ordering. More precisely, the staggered antiferromagnetic correlation develops along 
the chains in the xy planes, however the coupling between the neighboring chains 
remains frustrated: Hence, the system cannot develop three-dimensional magnetic 
long-range order. We examine this mean-field picture by using Monte Carlo simu- 
lation. 

Figure ITU1 shows the temperature dependence of the specific heat per site, which 
is calculated by fluctuations of the internal energy as 



The brackets denote thermal averages by Monte Carlo sampling. The specific heat 
shows a jump at To — 0.23 J, indicating a first-order transition. 

At this transition temperature, orbital ordering takes place. Figure ^2 shows 
the sublattice orbital moment, which is calculated in the form 
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Fig. 10. Temperature dependence of the specific heat per site [Eq. 14-71 1 for the model with nearest- 
neighbor interactions only. 



where the summation is taken over the sites within one of the four sublattices in 
tetrahedron unit cell. Here, Ii is the three-state clock vector at the site i which 
describes three different orbital states as shown in the inset of Fig, lilt Jj = (1, 0) for 
(dy Z ,d xy ), Ii = (-1/2, y/3/2) for (d yz ,d zx ), and Ii = (-1/2,-^/3/2) for (d zx ,d xy ) 
orbital occupations, respectively. It is found that the values of Mo for four different 
sublattices have the same value within the error bars so that we omit the sublattice 
index in Eq. (|4-8[) . As shown in Fig. 1111 Mo shows a clear jump at To — 0.23 J 
where the specific heat jumps. This indicates that a four-sublattice orbital ordering 
occurs at To- 

Figure WI\ plots the orbital distribution for four sublattices in tetrahedron unit 
cell, which is defined as 

n a = 1 r^- V (n ia ), (4-9) 

iVgjte * — * 

igsublattice 

where a = 1 (d yz ), 2 (d zx ), and 3 (d xy ). The results indicate that the orbital 
distributions suddenly change from equally distributed state n a ~ 2/3 in the high- 
temperature para phase to almost polarized state n a ~ or 1 below To- In the 
orbital ordered phase, d yz (a = 1) and d xy (a = 3) orbitals are occupied in the 
sublattices 1 and 4, and d zx (a = 2) and d xy (a = 3) orbitals are occupied in the 
sublattices 2 and 3. This orbital ordering structure is consistent with the mean-field 
prediction shown in Fig. EJa). 

Accompanying the transition at To, a tetragonal lattice distortion occurs dis- 
continuously. In Fig. 1131 we plot the average of the Jahn-Teller distortions which is 
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Fig. 11. The sublattice orbital moment in Eq. 14-811 . The inset shows the three-state clock vector 
for the orbital state. 



1.0 

0.8 

0.6 

0.4 

0.2 

0.0 
1.0 



(a) sublattice 1 



x n 3 




(b) sublattice 2 



(c) sublattice 3 



0.8 
0.6 - 
0.4 - 
0.2 

0.0 I AA/wwW ^ ^AAAA M fli iaa: ^ I 

0.0 0.1 0.2 0.30.0 



- (d) sublattice 4 





0.3 



Fig. 12. Electron density in each orbital for four sublattices 1-4 defined in Eq. 14-91 1. The data are 
calculated at L — 8. 
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The positive value of Q below To corresponds to a ferro-type tetragonal distortion 
with a compression of VOq octahedra. 
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Fig. 13. The average of the Jahn-Teller distortion defined in Eq. 14-1011 . 



Therefore, the discontinuous transition at To — 0.23J is ascribed to the orbital 
ordering with the tetragonal Jahn-Teller distortion. The orbital and lattice ordered 
state in the low-temperature phase is consistent with the mean-field prediction. 

So much for the orbital and lattice state — what about the spin sector? Monte 
Carlo results clearly indicate that spin frustration remains even below To as pre- 
dicted in the mean-field picture. Figure ITU plots temperature dependences of spin 
correlations for nearest-neighbor and third-neighbor pairs. The spin correlations 
are measured along the xy, yz and zx chains, respectively, as 



s {v) 



E 



(Si ■ Sj 



S. 



M 

3rd 



E 

«y»ei> 



(Si ■ Sj 



(4-11) 



where v = xy, yz, zx, and is the number of the bonds in the summation. The 
summations with (ij) and ((ij)) are for the nearest-neighbor sites and the third- 
neighbor sites along the v chains, respectively. 

As shown in Fig. 114( a). along the xy chains, the antiferromagnetic correlations 
develop below To- On the other hand, the spin correlations are strongly suppressed 
in the yz (and zx) directions as shown in Fig. 114( b). These results support the 
mean-field picture, which predicts the antiferromagnetic correlations along the xy 
chains and strong frustration between the chains. Since one-dimensional Heisen- 
berg spin chain cannot order at any temperature because of strong fluctuations, the 
model (|4-6|) . which includes only the nearest-neighbor couplings, does not show any 
magnetic ordering at finite temperature. This is the reason why we do not see any 
anomaly associated with magnetic ordering in the specific heat in Fig. I1UI 

4.2.3. Effect of third-neighbor superexchange interaction 

The mean-field analysis in Sec. 14.11 predicts that the magnetic transition in 
experiments is reproduced by including third-neighbor superexchange interactions. 
Here, we examine this prediction by using Monte Carlo technique. Now the total 
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Fig. 14. Spin correlations for the nearest-neighbor sites and the third-neighbor sites defined by 
Eq. I|4-ll[l along the (a) xy chains and (b) yz (or, equivalently zx) chains. The data are calculated 
for the model with nearest-neighbor interactions only with the system size L = 8. 



Hamiltonian reads 



77 77-ui] , rj3rd 
H - M SE + H SE + -".FT- 



Here the third-neighbor superexchange term is given by 



£J~3rd 



J 



«<i» 



o— AF + "o-F J ' 



(4-12) 



(4-13) 



where the summation is taken over the third-neighbor sites along the chains, and 

%-AF(F) are gi ven by Eqs. (j3-5f) and (|3-6|l . 

When we switch on the third-neighbor interaction J3, additional anomaly ap- 
pears in the specific heat except for the jump associated with the orbital and lattice 
orderings. Figure IT51 shows the specific heat per site at J 3 = 0.02 J. In addition to 
the jump at Tq ~ 0.19 J, a peak appears at ~ 0.115 J. This peak, not a jump, in- 
creases its height as the system size L, suggesting a divergence in the thermodynamic 
limit — a characteristic feature of a second-order phase transition. 

The jump at To — 0.19 J is ascribed to the discontinuous transition of orbital and 
lattice orderings. The order parameters behave in a similar manner as in Figs. ^2 
and El although the transition temperature Tq slightly decreases. On the other 
hand, the anomaly at Tn — 0.1 15 J is ascribed to a magnetic ordering as explained 
in the following. 

Figure [TdI shows the spin correlations in Eq. (|4-11|) for the model (|4-12|) at J3 = 
0.02 J. Compared to the results at J3 = in Fig. 1141 a big difference is found in the 
third-neighbor spin correlation along the yz (and zx) chains; the antiferromagnetic 
correlation rapidly develops below Xn ~ 0.1 15 J. This suggests an emergence of 
three-dimensional magnetic ordering there. 
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Fig. 15. Temperature dependence of the specific heat per site [Eq. l|J2j] for the model with the 
third-neighbor superexchange interaction J3 = 0.02J. 
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Fig. 16. Spin correlations [Eq. 14-lljl ] for J3 = 0.02J. The data are calculated at L — 12. 
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We calculate the order parameters for the magnetic ordering. As described in 
Sec. 14.11 in the mean-field level, three-dimensional magnetic ordering takes place on 
two independent sublattices which consists of (110) chains and (110) chains, respec- 
tively [See Fig. Efa)]. Here, we measure the staggered moment of these sublattice 
magnetizations, which is formally defined by 



M s 



(4-14) 
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Fig. 17. (a) The staggered moment defined in Eq. (b) The staggered magnetic susceptibility 

measured by fluctuations of Ms. The results are calculated at J3 = 0.02 J. 

where the form factor gi is given by 

gi = cos[27r(xj + yi)} + icos[27r(a?i - yi)}. (4-15) 

Here X{ and yi are the x and y coordinates of the site i measured in the cubic unit, 
respectively. Note that gi is specified only by the x and y coordinates of the site i 
because the z coordinate is uniquely determined within the cubic unit cell due to 
the special structure of the pyrochlore lattice. 

Figure lTTl plots the temperature dependences of the sublattice magnetization Mg 
in Eq. (|4-14|) and the corresponding susceptibility xs- As shown in Fig. 117( a). the 
staggered magnetization develops continuously below Tn ^ 0.115J and approaches 
the saturation moment Ms = 1 as T —* 0. The susceptibility shows a diverging 
behavior at Tn as shown in Fig. 117( b) . These results indicate that the phase transi- 
tion at Tn is a second-order one with the magnetic ordering of staggered sublattice 
moments. 

In the mean-field level, two sublattice moments are independent. It is suggested 
that quantum fluctuations can align two moments and select a collinear magnetic 
ordering as discussed in Sec. 14.11 While the effect of quantum fluctuations is beyond 
our classical approach, our Monte Carlo study retains thermal fluctuations that are 
completely neglected in the mean-field argument, and it is known that in many frus- 
trated systems thermal fluctuations also favor a collinear state. 15 - 1 Here we measure 
the collinearity between two sublattice moments Mi and M2 by 

2 L\(Mi) 2 (M 2 ) 2 / 3J 2 V 111 3 J' V 7 

where 612 is the angle between M\ and Mi- Note that P becomes 1 when the AF 
order is collinear, i.e., M\ || M2, and that P becomes when M\ and M2 are 
independent. 

Figure ITH1 plots the measure of collinearity P. Below Tn — 0.115 J, P becomes 
finite and grows rapidly; it approaches 1 as the system size increases. This indicates 
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Fig. 18. The measure of collinearity defined in Eq. (14-1611 at J3 = 0.02J. The lines are guides for 
the eyes. 



( d zx, d xy) 




Fig. 19. Monte Carlo results for (a) orbital-ordering pattern and (b) spin-ordering pattern in the 
lowest-temperature phase. In (a), filled (open) circles denote the sites where d zx and d xy (d yz 
and d xy ) orbitals are occupied. 

that the magnetic state is collinear below Tn- Consequently, the magnetic order is 
indeed the same one observed in experiments, that is, up-down-up-down- staggered 
pattern along the xy chains and up-up-down-down- four-times period along the yz 
and zx chains — see Fig. I19f b), 

Therefore, the present model (|4-12|) . which includes the third-neighbor superex- 
change interactions, successfully explains the experimental results for both struc- 
tural transition and antiferromagnetic transition. The structural transition at ~ 
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Fig. 20. Phase diagram for the spin-orbital-lattice coupled model determined by classical Monte 
Carlo calculations. The orbital and lattice transition at To (the dashed line) is a first-order 
transition, while the antiferromagnetic transition at Tn (the solid line) is a second-order one. 
The lines are guides for the eyes. 



50K is understood as the orbital ordering accompanied by the tetragonal Jahn- 
Teller distortion of a compression of VOg octahedra. The orbital ordering pattern 
is a layered type — the so-called A-type order, which is alternative stacking of 
(dyz, d X j / )-occupied layers and (d zx , d^-occupied layers in the z direction, as shown 
in Fig. ll9f a). The transition at ~ 40K is the antiferromagnetic transition. The tran- 
sition temperature Tn is lower than To because the spin frustration is reduced Below 
Tn, the complicated antiferromagnetic ordering appears as shown in Fig. ll9f b ), which 
is consistent with the neutron scattering result in Fig. |2J This magnetic ordering is 
stabilized by the third-neighbor superexchange J3 as well as thermal fluctuations. 
The latter fluctuation effect is important to realize the collinear spin configurations 
by order-by-disorder mechanism. 

4.2.4. phase diagram 

The antiferromagnetic ordering is stabilized by J3, and hence, the transition 
temperature Tn increases as J3. On the other hand, the orbital and lattice transition 
temperature To slightly decreases as J3 increases as found in previous sections. What 
happens for larger J3? 

We determined the phase diagram by changing J3. The result is shown in Fig. l2l)l 
The transition temperatures are determined by the specific heat, and the orbital, 
lattice and magnetic order parameters. As shown in Fig. I20[ for large values of J3 > 
0.03 J, Tn coincides with To, and there is only one discontinuous transition where 
both orbital and antiferromagnetic spin moments becomes finite discontinuously. 
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In the region of small J3 , our model shows two transitions at different tempera- 
tures as observed in vanadium spinels. By using estimates of parameters J ~ 200K 
and J3 ~ 0.02 J ~ 4K, 2 ) the transition temperatures in our results are To ~ 40K and 
Tn ~ 20K. The experimental values are ~ 50K and ~ 40K, respectively, in ZnV204. 
The semiquantitative agreement of the transition temperatures is satisfactory. 

4.3. Effects of quantum fluctuations 

We have determined the ordering pattern of spin configuration based on the spin- 
orbital-lattice model and shown that it is consistent with the experimental result of 
elastic neutron scattering. The size of the static magnetic moment, however, needs 
further consideration, since it is subject to fluctuations and renormalized to a smaller 
value. To compare with the value at low temperatures in experiment, we now show 
that quantum fluctuations are very large in frustrated pyrochlore lattice, particularly 
when the third-neighbor exchange coupling J3 is small, and the reduction of the 
moment is prominent. 

Let us consider the Heisenberg model, 

#spin = JijSi ■ Sj + J3 Si ■ Sj, (4" 17) 

(ij) ««» 

and calculate the moment reduction at zero temperature by using the linearized 
spin-wave theory. Here, the first sum is taken over nearest-neighbor pairs, and the 
exchange constant is Jij = J xy or Jjj = J yz ,zx depending on the bond direction due 
to the underlying orbital ordering in Fig. I19f a). J3 is the third-neighbor exchange 
coupling, shown in Fig. |Hfb). Note that J xy 3> J3 are both antiferromagnetic, while 
Jyz,zx is ferromagnetic and \J yz ,zx\ <C Jxy Our calculations follow a standard pro- 
cedure and the only special point in the calculations is a frustrated geometry of 
pyrochlore lattice and the anisotropic pattern of exchange couplings. The starting 
classical state is the one determined in the previous section and shown in Fig. 119( b). 
The magnetic unit cell has eight spins and we choose a pair of neighboring two 
tetrahedra (See Fig. I21|) . Quantum spin fluctuations around classical values are 
represented by the Holstein-Primakoff transformation in terms of boson operators. 
After inserting this for the Heisenberg model (|4T7|) . we keep terms up to bilinear in 
boson operators and neglect higher-order terms. This bilinear boson Hamiltonian 
contains terms of boson-pair creation and annihilation and it is diagonalized by a 
Bogoliubov transformation. New boson operators defined in the Bogoliubov trans- 
formation describe magnetic excitations from the renormalized ground state. The 
moment reduction is represented by energy dispersion of these elementary excita- 
tions and matrix elements with original bosons introduced in the Holstein-Primakoff 
transformation, and is calculated in this way. 

We show the moment reduction AS in Fig.^Uas a function of the third-neighbor 
exchange coupling J3. Here, we use J xy = 1 along the xy chains and J yz , zx = —0.1 
along the yz or zx chains. The renormalized size of static magnetic moment is then 
given as Mg = g^iS — AS), where g and /ib are g-factor and Bohr magneton, 
respectively, and we set an ordinary value g = 2, since the orbital magnetic moment 
is now quenched. As J3 decreases, AS grows and diverges slowly in the J3 — > limit. 
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Fig. 21. Magnetic unit cell in the spin ordered phase. 
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Fig. 22. Moment reduction as a function of J3. 



We show that this divergence is logarithmic in the following. 

An essential point is the presence of zero modes at ^3=0. We show the dis- 
persion relation of the lowest magnon mode along symmetric axes in the Brillouin 
zone at J3 = and 0.02 in Fig. 12^1 The magnon energy vanishes in two planes of 
k x = k y and k x = —k y as indicated in Fig. l2*5f a'). Therefore, the magnon dispersion 
is one-dimensional at J3 = in this sense. As we discussed in a previous section, 
chains in the xy planes have large antiferromagnetic interchain couplings (J X y) and 
these one-dimensional chains are weakly coupled with other chains in adjacent xy 
planes. Despite the presence of the weak ferromagnetic interplane couplings (J yz ,zx), 
magnons do not propagate along the z direction when k x = ±k y . Each spin in xy 
chain is connected to two spins in a chain in an adjacent xy plane. Quantum inter- 
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Fig. 23. The lowest magnon mode at (a) J3 = and (b) J3 = 0.02. 



ference of the two neighbor spins completely suppresses the propagation of magnon 
at these wavevectors, and this is also a manifestation of geometrical frustration ef- 
fects of pyrochlore lattice. In the vicinity of k x = ±k y planes, the magnons of the 
lowest mode have a linear dispersion of energy, cok oc \Sk\, where \5k\ is the dis- 
tance from the plane of k x = ±k y . This leads to a logarithmic divergence in the 
moment reduction, AS oc log (System Size). When the third-neighbor coupling J3 
is switched on, the zero modes acquire a finite energy that is proportional to J3 as 
shown in Fig. !23f b N ). Therefore, the logarithmic divergence is cut off at this energy 
scale, resulting in AS ~ log(l/J3). From Fig. for realistic value of J3 = 0.02 J, 
the moment reduction becomes AS ~ 0.5. This leads to a largely reduced staggered 
moment Ms ~ 1[/^b]) which is consistent with experimental results. 19 ** 



§5. Discussion 



5.1. Effect of relativistic spin-orbit coupling 

Thus far, we neglected relativistic spin-orbit coupling, and employed the model 
that describes an interplay between the spin-orbital superexchange interactions and 
the cooperative Jahn-Teller coupling of tetragonal symmetry. We showed that the 
model successfully explains many experimental results of thermodynamic properties 
in vanadium spinel oxides. After our works, there have been complemental studies 
which theoretically examine effects of the relativistic spin-orbit coupling. Here, we 
briefly discuss about them in comparison with our work. 

Tchernyshyov proposed a model assuming the dominant role of the relativistic 
spin-orbit coupling. 16 ) In his model, an effective angular momentum L' with the 
magnitude L' = 1 is introduced to describe (£23) 2 configurations. The orbital angular 
momentum is given by L = —L', and couples to the 5 = 1 spin S at each site 
antiferromagnetically, and consequently, the total momentum J' = L' + S with 
the magnitude J 1 = 2 is a conserved, good quantum number in the system. The 
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lowest energy state is the J' = 2 quintuplet. When the Jahn-Teller effect is included 
as a perturbation, a degeneracy remains partially in this manifold: The ground 
state energy is the same between J' z = state with an elongation of tetrahedron and 
J' z = ±2 states with a compression of tetrahedron. Moreover, this degeneracy cannot 
be lifted even when the superexchange interactions are included. Consequently, 
the model fails to reproduce by itself the spin and orbital ordered state observed 
in experiments. By putting as an input the experimental fact that the magnetic 
ordering takes place, the model predicts the ground state with J' z = ±2 and a 
compression of tetrahedron. 

The orbital ordering predicted in this way is different from our result. In Tch- 
ernyshyov's model, the orbital occupation becomes uniform in space, that is, at all 
sites 



The orbital state is a complex one, \d yz ± id zx ) at all sites. On the other hand, our 
model predicts the A-type orbital ordering with a staggered occupation of d yz and 
d zx orbitals in the z direction as shown in Fig. I19f a). As pointed out in Ref. 16), the 
crystal symmetry is different between these two orbital states. The A-type ordering 
in Fig. I19f a) breaks mirror reflections in (110) and (110) planes as well as diamond 
glides in (100) and (010) planes, on the other hand, the uniform ordering in Eq. 1)5-1(1 
preserves these symmetries. 

Experimentally, the lattice symmetry below the structural transition tempera- 
ture has not been completely determined yet. X-ray scattering results for polycrystal 
samples suggest the persistence of the mirror and glide symmetry at the low tem- 
perature phase. 17 )' 18 ) On the contrary, recent synchrotron X-ray data for a single 
crystal sample suggest a breaking of the glide symmetry although they cannot de- 
cide whether the mirror symmetry is broken or not. 19 ) This controversy probably 
comes from a small electron-phonon coupling in this ti g system. In order to settle 
the theoretical controversy, sophisticated experiments with large single crystals are 
highly desired. 

Another model has been proposed by Di Matteo et al. 13 ^ The model takes into 
account the spin-orbital superexchange interaction and the relativistic spin-orbit 
coupling on an equal footing. Mean-field analysis was applied to obtain the ground 
state. In the absence of the relativistic spin-orbit coupling, the ground state is 
the one shown in Fig. Efb) for finite Hund's-rule coupling r\ > as discussed in 
Sec. 13.21 When the relativistic spin-orbit coupling is turned on, the orbital state 
changes into a highly antisymmetric one (four sites in tetrahedron unit cell take 
different orbital states including both real and complex ones) with keeping the 3- 
up and 1-down spin configuration. Finally, for moderate values of the relativistic 
spin-orbit coupling, the ground state shows a transition into the spin and orbital 
ordered state which is compatible with Tchernyshyov's prediction by assuming the 
dominant relativistic spin-orbit coupling. In this sense, this study may interpolate 
our model and the Tchernyshyov's model. In their treatment, however, the Jahn- 
Teller coupling is completely neglected. As we have seen in Sec. 14.11 the mean- 
field ground states are drastically changed by the effects of the Jahn-Teller lattice 
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distortion. It is interesting to investigate with including the cooperative aspect of 
the lattice distortions. 

5.2. Direction of magnetic moment 

Experimentally, the spontaneous magnetic moment appears along the z direc- 
tion. This is naturally explained by the scenario based on the relativistic spin-orbit 
coupling in the previous section as follows. The complex orbital state \d yz ± id zx ) 
has a finite matrix element of the z component of orbital angular momentum L. 
This aligns the magnetic moment along the z direction through the spin-orbit cou- 
pling. On the other hand, the spin part of our model is completely isotropic, and 
no direction is preferred for spin ordering. Therefore, one needs to introduce a pro- 
cess with spin anisotropy to explain the observed spin direction. Here, we discuss 
that the direction of moment is reproduced by our model when we include quantum 
fluctuations in the orbital sector as well as the relativistic spin-orbit coupling in a 
perturbative way. 

As described in Sec. 13.11 the orbital part of our model is completely classical, 
three-state clock type interaction, because we take into account the dominant a- 
bond transfers only. If we include small contributions from tt and <5-type transfer 
integrals, the orbital interaction is no longer classical, that is, contains small off- 
diagonal exchange interactions. This should lead to quantum fluctuations in the 
orbital sector. Our results in Sec. ^ showed that d xy orbital is lowered and occupied 
at every site, and that d yz and d zx orbitals are alternatively occupied as long as the a 
contribution is dominant. Under this situation, we expect that the orbital quantum 
nature, coming from small tt and 5 contributions, leads to fluctuations between d yz 
and d zx orbitals dominantly. This d yz —d zx fluctuation gives rise to a finite matrix 
element of L z . As mentioned above, L couples to the spin S antiferromagnetically 
once we consider a finite relativistic spin-orbit coupling, and therefore, the induced 
L z component aligns the magnetic moment in the z direction. 

5.3. Weakly- coupled S = 1 chains — Haldane gap? 

As discussed in Sec. @J the orbital ordering reduces the magnetic frustration in 
our model. Below the orbital ordering temperature, antiferromagnetic correlations 
develop along the xy chains, and the coupling is frustrated between neighboring 
chains. Hence, the system can be regarded as weakly-coupled S = 1 spin chains 
under orbital ordering. It is well known that the ground state of antiferromag- 
netic S = 1 chain is spin-singlet, and a gap opens between the ground state and 
the lowest triplet excitation — the so-called Haldane gap. 20 ) This is purely quan- 
tum mechanical effect, and hence it is beyond our present theory by mean-field 
argument and Monte Carlo simulation to examine the possibility of Haldane gap. 
Experimentally, vanadium oxides exhibit magnetically-ordered ground state, not the 
spin-singlet state, and we here address this issue in the light of our results and make 
some remarks. 

It is known that the Haldane gap is fragile against the interchain couplings. For 
example, when antiferromagnetic 5 = 1 chains are connected to form of a simple 
cubic lattice structure, the critical value of the interchain exchange coupling is esti- 
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mated at 2 — 3% of the intrachain exchange coupling. 21 - 1 In our pyrochlore system, 
the relevant interchain coupling is the third-neighbor exchange J3 in Fig. |HJa). The 
coupling J3 constitutes a complicated three-dimensional network without any frus- 
tration among the chains. In addition, the neighboring xy chains are coupled by 
weak ferromagnetic coupling as well as second-neighbor coupling although they are 
frustrated due to the pyrochlore geometry. We expect for realistic parameter regime 
that the Haldane gap state is unstable against these interchain couplings. 

On the other hand, as shown in the phase diagram in Fig. I2U1 the transition 
temperature for the magnetic ordering decreases rapidly as J3, and becomes zero in 
the limit of J3 within our classical model. This implies that once all farther-neighbor 
couplings are irrelevant, a competition occurs between the long-range magnetic order 
and the quantum Haldane gap formation, leading to a quantum critical point. This 
interesting problem is left for further study. 

§6. Summary 

In this paper we have reviewed our recent studies on the geometrically frus- 
trated systems, vanadium spinel oxides ^4V204 with divalent ^4-site cations. These 
are spin-1 antiferromagnets on the frustrated pyrochlore network of vanadium sites. 
We have employed an exchange model of electron spin and orbital degrees of freedom 
with including coupling to Jahn-Teller lattice distortion. By using mean-field ap- 
proach and numerical calculations, we have shown that this model can describe two 
phase transitions in this system. The structural transition at a higher temperature 
is ascribed to orbital ordering accompanied by lattice distortion. The orbital order- 
ing leads to an spatially anisotropic pattern of spin exchange couplings. Geometrical 
frustration of spins is thus relaxed, and this induces a magnetic transition. Spin con- 
figuration in the magnetic phase was determined and we found that it agrees with the 
one observed by neutron scattering experiments. It was also found that the third- 
neighbor spin exchange couplings are important to determine the spin configuration 
and the magnetic transition temperature. A realistic value of the third-neighbor 
coupling gives a size of the magnetic moment and a transition temperature, both of 
which are quite well in agreement with the experimental values. These semiquan- 
titatively good agreements with experimental results support our choice of model 
Hamiltonian and also our picture of the two transitions. A few points, however, 
remain to be investigated in more details. One is more systematic analysis of the 
Jahn-Teller effects. So far, we have considered only uniform modes of lattice dis- 
tortion, experimentally observed ones, but it is important to check that the uniform 
mode is the most unstable one. Another point is more quantitative analysis of the 
effects of transverse orbital fluctuations. These are introduced by electron hoppings 
in 7r- and <5-bonds. The orbital part of the effective exchange model remains strongly 
anisotropic of three-state Potts type even with including transverse orbital couplings, 
but it is interesting to investigate the effects of these couplings. The transverse cou- 
plings induce orbital angular momentum, and therefore the relativistic spin-orbit 
coupling may show interesting effects. 
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